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ABSTRACT 

We present a model for microstructure in pulsar radio emission. We propose that 
micropulses result from the alteration of the radio wave generation region by nearly 
transverse drift waves propagating across the pulsar magnetic field and encircling 
the bundle of the open magnetic field lines. It is demonstrated that such waves can 
modify significantly curvature of these dipolar field lines. This in turn affects strongly 
fulfillment of the resonance conditions necessary for the excitation of radio waves. The 
time-scale of micropulses is therefore determined by the wavelength of drift waves. 
Main features of the microstructure are naturally explained in the frame of this model. 
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1 INTRODUCTION 

Intensity variations of pulsar radio emission have several dif- 
ferent time-scales, starting from a few nanoseconds (nanos- 
tructure). Single pulses consist of subpulses (there are usu- 
ally about 3 -f 5 of them in a single pulse window) which of- 
ten undergo drift, i.e. transition in phase. The nature of sub- 
pulses and most of their various features (such as circular po- 
larisation, drifting subpulses, mode switching, nuUings and 
phase memory phenomenon) were explained by the Plasma 
model for pulsar emission developed in the series of papers 
(see, e.g. Kazbegi et al. 1987, 1991a, 1991b, 1991c, 1996). 
In this paper we present an explanation of the pulsar mi- 
crostructure in the frame of the Plasma model. 

Let us first summarise the main observational features 
of microstructure. Micropulses are ultrashort intensity vari- 
ations within individual pulses with the following properties. 
(1) Characteristic time-scale is 20 30 ^s^ with an upper 
limit of the time width distribution 0.1 -i- 1 ms (Boriakoff 
1996). (2) Individual micropulse widths are constant with 
frequency. (3) Micropulse widths are about twice smaller 
than the distance between micropulses. (4) The time sepa- 
ration (phase) of any micropulse with respect to the fiducial 
point is constant with frequency, i.e. it is simultaneous at all 
frequencies at source (let us note here that the position of 
subpulses typically varies with the frequency). (5) In about 
25 per cent of the cases micropulses are observed in quasi- 
periodic sequences (trains) which can be detected simulta- 
neously at widely separated frequencies, e.g. at 430 MHz 
and 1.4 GHz. (6) The typical life-time of a micropulse does 



Time-scale generally ranges from 1 to 1 ms (Cordes 1979). 



not exceed one second. (7) Modulation depth depends on 
the frequency as u~^'^. 

Lange et al. (1998) reported on observations of mi- 
crostructure at 4.85 and 1.41 GHz for a few pulsars. It ap- 
peared that for all observed pulsars a large fraction of the 
single pulses (varying between 30 and 70 per cent) show 
microstructure. For the pulsars where also low frequency re- 
sults are available (Popov et al. 1988), there is no significant 
diflterence between microstructure properties at high and low 
frequencies. These results confirmed that microstructure is 
a common property of the pulsar emission and not only an 
additional feature of a few strong pulsars. Thus, the mi- 
crostructure represents perhaps one of the fundamental fea- 
tures of pulsar pulses. Each self-consistent theory of pulsar 
radio emission should be able to explain existence of this 
phenomenon. 

Two general types of models were presented for pul- 
sar microstructure, often regarded as beaming models and 
temporal models, respectively (Chian & Kennel 1983). In 
beaming models the observer's line of sight sweeps across a 
non-uniform pulsar beam, which results in rapid intensity 
fluctuation of observed radiation (Benford 1977; Ferguson 
1981). In the angular beaming model proposed by Cordes 
(1979) and Gil (1982, 1986) both subpulses and micropulses 
are generated by the same narrow-band emission mecha- 
nism. Here micropulses correspond to thin plasma columns 
flowing along dipolar magnetic field lines, and the width 
of micropulses is determined by relativistic beaming. The 
model explains frequency stability of micropulses, in con- 
trary to subpulses, basing on the geometrical consideration. 
The temporal models assume that the pulsar radiation is 
modulated while propagating through the magnetospheric 
plasma (see, e.g. Harding & Tademaru 1981; Chian & Ken- 
nel 1983). 



© 1999 RAS 



2 Machabeli et al. 



Basing on the similarity of the microstructure time- 
scales over the broad frequency range, Lange et al. (1998) 
claimed that the micropulse duration is generated by the di- 
mension of the emitting structure. Below we argue that the 
microstructure is caused by alteration of radio wave gen- 
eration region by nearly transverse drift waves propagat- 
ing across the magnetic field and encyrcling the open field 
lines region of the pulsar magnetosphere. Our mechanism 
explains naturally important features of microstructure. 



2 EMISSION MODEL 

It is generally assumed that the pulsar magnetosphere is 
fiUed by dense relativistic electron-positron plasma flow- 
ing along the open magnetic field lines, which is gener- 
ated as a consequence of the avalanche process first de- 
scribed by Goldreich & Julian (1969) and developed by 
Sturrock (1971). This plasma is multicomponent, with a 
one-dimensional distribution function (see Fig. 1 in Arons 
1981), containing: (i) electrons and positrons of the bulk 
of plasma with mean Lorentz factor of 7p and density n^; 
(ii) particles of the high-energy 'tail' of the distribution 
function with 74 and nt, stretched in the direction of pos- 
itive momenta; (iii) the ultrarelativistic (7b ~ 10®) pri- 
mary beam with so called 'Goldreich- Julian' density nt ~ 
7 X W~'^BoP^^{Ro/r)^ [cm"^] (where P is a pulsar period, 
Ro is a neutron star radius, Bq is a magnetic field value at 
the stellar surface and r is a distance from the neutron star's 
centre), which is much less than Up {k = Up/ub ~ 10^~®). 
Such a distribution function should generate various wave- 
modes in certain conditions. These waves then propagate in 
the pair plasma of a pulsar magnetosphere, transform into 
the vacuum electromagnetic waves as the plasma density 
drops, enter the interstellar medium, and reach an observer 
as the pulsar radio emission. 

An important feature of the pulsar radio emission is 
that the circular polarisation is observed in the vicinity of 
the maximum of a pulse (see, e.g. Taylor & Stinebring 1986). 
At the same time, the exact maximum of the pulse power 
corresponds to the zero circular polarisation (as the latter 
changes its sign) . This means that the waves leave the mag- 
netosphere propagating at relatively small angles to the pul- 
sar magnetic field (Kazbegi et al. 1991b). 

An extensive analysis have been conducted (Volokitin 
et al. 1985; Arons & Barnard 1986; Lominadze et al. 1986) in 
order to study the dispersion properties of the waves prop- 
agating through the highly magnetised relativistic electron- 
positron plasma of pulsar magnetospheres. In the general 
case of oblique propagation with respect to the magnetic 
field three different wavemodes can be distinguished, these 
are: (i) the purely electromagnetic X-mode (also called t- 
wave), (ii) the subluminous Alfven mode (it-wave) and iii) 
the superluminous O-mode (L-wave). The last two modes 
are of mixed electrostatic-electromagnetic nature. Electric 
field vectors B° and of the O and A-modes lie in the 
(kB) plane, while the electric field of the X-mode E'^ is 
directed perpendicularly to this plane. Here k a wave- vector 
and B is a local magnetic field. 

Particles moving along the curved magnetic field un- 
dergo drift transversely to the plane of field curvature, with 
the velocity 



(1) 



where Ug = (eB/mc), Rc is the curvature radius of the 
dipolar magnetic field line, 7 is the relativistic Lorentz factor 
of a particle, and v^p is a particle velocity along the magneic 
field line. Here and below the cylindrical coordinate system 
{x, r, (fi) is chosed, with the x-axis directed transversely to 
the plane of the magnetic field curvature, while r and ip are 
the radial and azimuthal coordinates, respectively. 

Generation of X and A-modes is possible after satisfac- 
tion of one of the following resonance conditions: 

ui — k^pV^p — kxU — — —, Cyclotron instability, and/or (2) 

7rcs 



ui — kipV^ — kxU — 0, Cherenkov instability. 



(3) 



These conditions are very sensitive to the parameters of the 
magnetospheric plasma, particularly to the value of the drift 
velocity (equation 1), hence to the curvature of magnetic 
field lines (Kazbegi et al. 1996). 

Let us notice that the Plasma model is in accordance 
with observed systematic increase of component separation 
and profile widths with decreasing frequency, often called 
radius-to- frequency mapping (RFM). The RFM implies that 
the generation region is restricted to narrow range of alti- 
tudes. Besides, low-frequency waves are generated at rela- 
tively lower altitudes within this region. Namely, provided 
that magnetic field is dipolar far from the stellar surface, 
the wave frequency scales with the distance as 



(4) 



in the case of the cyclotron (equation 2) mechanism, and as 

'Ro' 



UJo tx 



(5) 



in the case of the Cherenkov (equation 3) mechanism 
(Kazbegi et al. 1991a). It is obvious that both mechanisms 
satisfy RFM, hence the emission region is localised to dis- 
tinct altitudes. 



3 GENERATION OF DRIFT WAVES 

It was shown (Kazbegi et al. 1991c, 1996) that in addition 
to the X and A-modes (whose characteristic frequencies fall 
into radio band) propagating with small angles to the mag- 
netic field lines, the very low-frequency nearly transverse 
drift waves can be excited. They propagate across the mag- 
netic field, so that the angle 6 between k and B is close to 
n/2. In other words, k_L/k^ S> 1, where k± = {k^ + k'^Y^'^ . 
Assuming 7(tj/c>j^) < 1, (ita/c)^ < 1, k^/k^ < 1 
and kr — > 0, dispersion equation of the drift wave writes 
(Kazbegi et al. 1991c, 1996) 



= i+y 5^ 



df 



kxUa 97 



d7, (6) 



where a denotes the sort of particles and Up^ = Annpae^ /m. 
Let us assume that 



to = kxUb + kuiVu, + a, 



(7) 
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where ut, is a drift velocity of the beam particles (equa- 
tion 1). Integration in parts and summation over the sorts 
of particles in equation (6), and use of equations (7) yields 



1 kxUp Ljf, kxUt kj^c 



(8) 



where subscripts 'p' and '&' denote values of the quanti- 
ties corresponding to the bulk of plasma and the beam, 
respectively. Note that a small term {loI /j^uja) has been ne- 
glected in equation (8). Thus, the frequency of a drift wave 
uo = Ke(jj writes 



LJo — kxUf, + k^v^ 



(9) 



where k^v^p <C k^ui,. 

Neglecting the second and the third terms in equa- 
tion (8), and solving this equation for the imaginary part 
of the complex frequency |a| ~ F = Imuj, we obtain that 
the growth rate of drift waves is a maximum when 



ki < 



In this approximation the growth rate writes 



1/2 



lb 



1/2 



kxUb- 



(10) 



(11) 



The drift waves propagate across the magnetic field and 
encircle the region of the open field lines of the pulsar magne- 
tosphere. They draw energy from the longitudinal motion of 
the beam particles, as in the case of the ordinary Cherenkov 
wave-particle interaction. However, they are excited only if 
kxUb 7^ 0, i.e. in the presence of the drift motion of the beam 
particles. Note that these low-frequency waves are nearly 
transverse, with the electric vector directed almost along 
the local magnetic field. 



4 THE MODEL FOR MICROSTRUCTURE 

Let us assume that a drift wave with the dispersion de- 
fined by equation (9) is excited at some place in the pulsar 
magnetosphere. It follows from the Maxwell equations that 
Br = Eip (kxc/ujo), hence B,- 3> E^p for such a wave. There- 
fore, excitation of a drift wave causes particular growth of 
the r-component of the local magnetic field. 

The field line curvature pc = is defined in a Carte- 

sian frame of coordinates as 



Pc 



1 + 



\dxj 



-3/2 2 
d y 



(12) 



where dy/dx — By/Bx- Using (V B) = and rewriting 
equation (12) in the cylindrical coordinates we get 



Pc = - 



IB^ 1 1 B; dBr 



r B r B B^ d<p 



(13) 



Here B = (B^ + B^f^^ ^ B^[l + {B^/2Bl)] . Assuming 
that fc^r S> 1 we obtain from equation (13) 



= r ( 1 - ^^^^ 



(14) 



From equation (14) it follows that even a small variation of 
Br causes significant change of pc (since k^r 3> 1). 



Thus, the drift wave affects significantly the curvature 
of magnetic field lines, which in turn alters the drift velocity 
(equation 1) of particles. On the other hand, the resonance 
conditions (2) and (3), which in the Plasma model are re- 
sponsible for the radio wave generation in pulsar magneto- 
spheres, are very sensitive to the parameters of the magneto- 
spheric plasma, particularly to the value of the drift velocity 
u. Therefore, any variation of the magnetic field line curva- 
ture affects strongly the process of waves excitation. It fol- 
lows that, within the generation region, the resonance con- 
ditions can be fulfilled only at the definite places ('emitting 
spots'), corresponding to the definite ('favorable') phases of 
the drift wave. The characteristic linear size Xm of a single 
emitting spot should be of order of the half-wavelength of 
the drift wave Xm ~ Xj_/2 — n/kj_. 

Assuming that the generation region corotates with a 
pulsar and introducing the altitude Rem of this region from 
the stellar surface, we can estimate the characteristic time 
during which the observer sweeps a single emitting spot 



An 



1 



2 k±Rc 



lO"'" {k^TZc 



(15) 



where Q = 2n/P is the angular velocity of a pulsar and 
Tlem = Rcm/RhC IS SL dimensionless altitude of the gener- 
ation region from the stellar surface measured in the units 
of the light cylinder radius 7?lc ~ 5 x lO^P [cm]. Let us 
calculate the drift wave-number in the generation region 
of a 'typical' pulsar (i.e., the pulsar with P = 0.5 s and 
P — 10~^^), for the magnetospheric parameters used in the 
Plasma model (7^ ~ 3 and k ~ 10^). Note that according 
to this model radio emission originates at TZcm ~ 0.5 ± 0.2. 
Thus, from equation (10) we obtain that 



fcx < 3 X 10 ^ cm"\ 



(16) 



^ cm ^, in agreement 



Speculating that kj_ 10"" ^ 10 
with equation (16), we obtain from equation (15) the char- 
acteristic times Tm ~ 10~^ -ir 10~* s, which agree well with 
the typical time-scales of the microstructure. Note that the 
case fcj_ ~ 10~^ cm~^ at TZcm ~ 0.5 corresponds to the 
typical observed micropulse width of about 20 /is, whereas 
drift waves with the longer wavelengths k± ~ 10~^ cm~^ 
account for the time-scales corresponding to the upper limit 
of the microstructure widths distribution. From the model 
presented above naturally follow all the other features of the 
microstructure. 

An observer distinguishes between the micropulses by 
their width which in turn is determined by the wave-length 
of the drift wave. On the other hand, the position of the 
emitting spots depends on the latitude, as it follows from 
equation (14). At the same time, the generation region, al- 
tered by the drift wave, is rather localised in both mecha- 
nisms (see equations 2 and 3), and the broad band of ra- 
dio waves with different frequencies is excited there. So the 
observer detects different frequencies which are coming al- 
most from the same place of the magnetosphere, i.e. from 
the same altitude and latitude. This explains the stability 
of the micropulse phase with respect to the pulse fiducial 
point at different frequencies. The quasi-periodic trains of 
micropulses observed at widely separated frequencies thus 
correspond to favorable phases of a drift wave. 

From equation (11) it follows that the growth rate of 
drift waves increases with the altitude from the stellar sur- 
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face. Hence, the alteration of the generation region by these 
waves is stronger at higher ahitudes where the lower fre- 
quencies are excited (equations 2 and 3). This explains why 
the probability of microstructure detection, as well as its 
modulation depth, decrease at higher frequencies. 

Let us notice that although <^ kx for the drift waves, 
there still exists some non-zero k^. From equation (9) it fol- 
lows that the phase velocity of the drift wave has a longi- 
tudinal component v^/c ~ k^p/kx which gradually draws 
it away from the resonance region. This implies that the 
corresponding micropulse should disappear after a period of 
time 



cm kx 
27T 



(17) 



where Ai?cm is a longitudinal size of the generation re- 
gion and ATZc-m is the same but in the units of the light 
cylinder radius. Assuming that A7?.cm ~ 0.4 we obtain for 
the 'typical' pulsar (with P = 0.5 s) that /Stm « 1 s if 
k^/kx ~ 0.03 ^ 1. This corresponds to the observed upper 
bound of the micropulse typical life-time. 



5 SUMMARY 

In this paper we attempt to explain the existence and various 
properties of the microstructure in the pulsar radio emission 
basing on the Plasma model. According to suggested mech- 
anism, the radio wave generation region is strongly altered 
by so called drift waves. These nearly transverse waves prop- 
agate across the local magnetic field, encyrcling the whole 
bundle of the open magnetic field lines, and cause significant 
change of field lines' curvature. This alters the curvature 
drift velocity (equation 1), therefore influencing fulfillment 
of the resonance conditions (equations 2 and 3) necessary 
for excitation of radio waves. Thus, radio waves are gener- 
ated only at the definite regions (emitting spots) placed at 
approximately the same altitude from the star but different 
latitudes, corresponding to definite phases of the drift wave. 
The time-scale of the micropulses is thus determined by the 
characteristic time needed to the observer's line of sight to 
sweep an emitting spot with characteristic transverse size 
being about half- wavelength of the drift wave (equation 15) . 
Characteristic times estimated from our model are in a good 
accordance with the typical observed time-scales of the mi- 
crostructure. 

Most of the important features of the microstructure are 
successfully explained by the model presented in this paper. 
However, we would like to stress that the results of the model 
are qualitative, and it does not pretend to strict quantitative 
explanation of all the existing observational data. 
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